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Abstract 

A new formalism is presented for analytically obtaining the probability den- 
sity function, Puis), for the distance between two random points in an n- 
dimensional sphere of radius R. Our formalism allows Pnis) to be calculated 
for a sphere having an arbitrary density distribution, and reproduces the 
well-known results for the case of a sphere with uniform density. The re- 
sults find applications in stochastic geometry, probability distribution theory, 
astrophysics, nuclear physics, and elementary particle physics. 
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I. INTRODUCTION 



In a recent paper geometric probability techniques were developed to calculate the 
probability density functions (PDFs) which describe the probability density of finding a 
distance s separating two points distributed in a uniform 3-dimensional sphere and in a 
uniform ellipsoid. Our focus in the present paper will be on the probability density functions 
Pn{s) for an n-dimensional sphere of radius R characterized by + + ■ ■ ■ + < R^, 
where xi, X2, and x„ are the corresponding Cartesian coordinates. (In the mathematical 
literature this is sometimes termed an n-dimensional ball). As discussed in Refs. |]T]-Q, 
these results are of interest both as pure mathematics and as tools in mathematical physics. 
Specifically, it was demonstrated in Ref. that geometric probability techniques greatly 
facilitate the calculation of the self-energies for spherical matter distributions arising from 
electromagnetic, gravitational, or weak interactions. The functional form of Pn{s) is known 
for a sphere of uniform density, and hence the object of the present paper is to generalize 
the results of Refs. 0-0 to the case of an arbitrarily non-uniform density distribution by 
using our method. As an application of these results we will consider the neutrino-exchange 
contribution to the self-energy of a neutron star modeled as series of concentric shells of 
different constant density. Other applications will also be discussed. 

In this paper we present a new technique for obtaining the analytical probability density 
function P„(s) for a sphere of n dimensions having an arbitrary density distribution. To 
illustrate this technique, we begin by deriving the PDF, Pn{s), for an n-dimensional uniform 
sphere of radius R, and compare our results to those obtained earlier by other means pj-Q . 
We then extend this technique to an n-dimensional sphere with a non-uniform but spherically 
symmetric density distribution. We explicitly evaluate the analytical probability density 
functions for certain specific density distributions, and then use numerical Monte Carlo 
simulations to verify the analytical results. 

Finally our formalism is generalized to an n-dimensional sphere with an arbitrary den- 
sity distribution, and leads to a general-purpose master formula. This formula allows one 
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to evaluate the PDF for a sphere in n dimensions with an arbitrary density distribution. 
After verifying that the master formula reproduces the results for uniform and spherically 
symmetric density distributions, we analytically evaluate the probability density functions, 
for 2, 3, and 4-dimensional spheres having non-uniform density distributions. The analytical 
results are then verified by the use of Monte Carlo simulations. The outline of this paper 
is as follows. In Sec. |T| we present our new formalism and illustrate it by rederiving the 
well-known results for a circle and for a sphere of uniform density. In Sec. |T| we extend 
this formalism to the case of non-uniform but spherically symmetric density distributions. 
In Sec. 1^ we develop the formalism for the most general case of an arbitrary non-uniform 
density distribution. In Sec. we present some applications of our formalism to physics. 
These include the mth moment (s™) for a sphere of uniform and Gaussian density distri- 
bution. Coulomb self-energy for a collection of charges, z/z/-exchange interactions, obtaining 
the probability density functions for multiple-shell density distributions found in neutron 
star models ^,0, and the evaluation of some geometric probabihty constants 

II. UNIFORM DENSITY DISTRIBUTIONS 

A. Theory for a uniform circle 

In this section we illustrate our formalism by deriving the PDF for a circle of radius 
R having a spatially uniform density distribution characterized by p = constant. For two 
points randomly sampled inside the circle located at ri and r2 measured from the center, 
define s = s{fi,f2) = r2 — ri and s = \s\. To simplify the discussion, we translate the 
center of the circle to the origin so that the equation for the circle is + = R^. It is 
sufficient to initially consider those vectors s which are aligned in the positive x direction, 
since rotational symmetry can eventually be used to extend our results to those vectors s 
with arbitrary orientations. We begin by identifying those pairs of points, ri and r2, which 
satisfy s = r2 — ri = sx, where < s < 2R. One set of points is obtained if the points 1 are 
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uniformly located on a line Li given by a; = — s/2 and the corresponding points 2 are on a 
line L2 given hj x = +s/2 as shown in Fig. ^ Notice that x = ±s/2 are two symmetric lines 
with respect to reflection about x = 0. Another set of points is obtained if n is uniformly 
located in the area Ai and f^2 is located in the area A2 as shown in Fig. 0(a). Note that 
Ai and A2 are congruent. The only remaining set of points that need to be considered are 
those for which fi is uniformly located in the area A3 while r2 is in the area A4 as shown 
in Fig. ^(b). As before ^3 and A^ are congruent. Furthermore, A4 is a reflection of A2 
about X = s/2, while A3 is a reflection of Ai about x = —s/2. As discussed in Appendix 0, 
the union of A2 and A4 {A2 U A4) is the overlap area between the original circle Co and an 
identical circle Ci whose center is shifted from (0,0) to (|s| ,0). Similarly, AiU^s is the 
overlap area between the original circle and an identical circle C2 whose center is shifted 
from (0,0) to (— |s|,0) as shown in Fig. ^. Since AilJ^s and A2U^4 are identical, it 
follows that the probability of flnding a given s = \sx\ in a uniform circle is proportional to 



A.\jA.=J,d.l^dy + l^d.l^__dy. (1) 



R rVR^-x'^ ^•a/R2-(x-s)2 

dx / dy = dx dy. 
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Notice that 



V/?2-a;2 Js-R J-^R2-(x-s)'^ 

Using rotational symmetry of Cq, this result can apply to any orientation of s between and 
27r, and hence the probability of finding a given s = \s\ is proportional to |s | Jq^ dcj) = 27rs. 
We denote this PDF for a uniform circle by P2{s), and impose the normalization requirement 

r2R 

/ P2is)ds = 1. (2) 
JO 

We then have 



27rs dx L^^^ dy + //_^^ dx dyj 

P2i-^) = r . (3) 

2.8 (^/f dx dy + dx /y^^=^ dyj ds 

Equation (^) can be simplified to 



lo^ {s jf dx dy) ds 

Equation (||) is identical to the results obtained in Refs. by other means. 

The conclusion that emerges from this formalism is that the probability of finding two 
random points separated by a vector s in a uniform circle can be derived by simply calculating 
the overlap region of that circle with an identical circle obtained by shifting the center from 
the origin to s. In the following sections, we show that this result generalizes to higher 
dimensions, and provides a simple way of calculating Pn{s) for n > 3. 

B. Theory for a uniform sphere 

The above formalism for a 2-dimensional uniform circle can be extended to a 3- 
dimensional uniform sphere with radius R. For a given s, we arbitrarily select the positive 
z direction and study the distribution of vectors s = r2 — ri in this direction. The areas Ai, 
A2, A3, and ^4 in Sec. [II A| are replaced by the volumes Vi, V2, V3, and V4. The PDF ^3(5) 



is therefore proportional to the overlap volume Vi U ^ = ^ U where 
V2 I J V4 = / dz dx / dy + dz dx / dy. 



(6) 



Notice that 



dz / dx / dy = / dz dx / dy. (7) 



Is-R J-^R^-{z-s)^ J-y/R2~{z~-s)^ 

In a 3-dimensional space, s can range from to vr in the 6 direction, and from to 2tt in 
the direction. For a given s with any orientation, the PDF ^3(5) is also proportional to 
Atts'^ by using rotational symmetry of a uniform sphere. Following the previous discussion, 
we thus arrive at the following expression for the PDF for a uniform sphere: 
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The result in Eq. (^) agrees exactly with the expression obtained previously in Refs. 



C. Representations , general properties, recursion relations, and generating 

functions of Pn{s) 

The present formalism can be readily generalized to obtain P„(s) for an n-dimensional 
sphere (n-sphere or n-ball) of radius R. From Eqs. and (|]) we find 



nl-^) = 7 T . 



(10) 



We can rewrite Eq. (|ToD as 



S" ^ //^ dXn C^a;! / ^ dx2 ■ ■ ■ / ^ /— ; — — dxn^i 



■2 — — /-2 , 



Pn{s) = , V, ^ '"^ (11) 



lo^ Is" 1 /f J^-ji dxi J ^ /— — ^ dx2--- !^ r-^ dxn-i 1 ds 







■^n-2 
( 






■^n-2 






■■■a;2 



where 



R=^R^- xl (12) 

As shown in Refs. P-pl]], the volume V{n, R) of an n-dimensional sphere of radius R is given 
by 

V (n, R) = dxi / , ax2 ■ ■ ■ / ax„ = — 7 ^. (13) 

^ ' ^ i-i? V-^]^ y_^fi2_,2_,2_...,2_^ " r(l + f) ^ ^ 

Equations. ( pUj ) and (|TT]) can then be reduced to a general simplified equation for an n- 
dimensional uniform sphere of radius R: 
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Qn{s) = I [r^-x^) ' dx, (22) 



2 

' dx. (23) 



It follows that the denominator of Eq. (14), which is the normalization constant, can be 
written as C{2R] 0,n) where 

r2R 



C{2R- 0,n) = / Tn{s)ds. 
Jo 



When n is an even integer, C{2R; 0,n) can be expressed in terms of the gamma and beta 
functions by noting that 

C{2R; 0, n) = il^ll!!^- ^ ALliMlll)^- = ^b(^ + \^-) i?-. (24) 
^ ' ' ^ 2n n!! 2n T (f + 1) 2n V2 2' 27 ^ ' 

Similarly, when n is an odd integer, 

C(2R; O.n) = i<^fi- = + U) (25) 



n nil 2n \2 2 2, 

We thus find that the normalization constant has the same functional form irrespective of 
whether n is even or odd, when expressed in terms of the beta function. 
If we introduce the variable t = R? — x^ and note that 

Q.{s) = -B^_^ + _B + _ _j + _j , (26) 

we can then rewrite Pni^s) in the form 

V2 ^ 2' 2^ 

where 



(p, q) is the incomplete beta function defined by 



8 



Bx (p, q) 



-1 



(30) 



and Ix is the normalized incomplete beta function defined by 

lAp. q)-^ r t (1 t) dt - ^^^^ . [31) 

Eq. ( [28| ) is the expression for P„(s) obtained in Refs. [0,^, and hence we have demonstrated 
that the classical results can be reproduced by using the formalism developed here. 

We find the following general properties of P„(s) and its derivative Pn{s) at the lower 
bound s = and at the upper bound s = 2R: Pi(0) = 1/R, P„(0) = forn > 2, P„(2P) = 0, 
P[{0) = -1/2R\ P^(0) = 2/P^ P;;(0) = forn > 3, Pi'(2P) = -1/2P^ and P;;(2P) = 
for n > 2. 

It is of interest to express the probability density functions Pn(s) in terms of generating 
functions from which several representations and recursion relations for P„ (s) can be derived. 
This can be achieved by observing that 



Pt 



n 1 1 



+ 



2' 2 



n\ \ dh 



F{h = 0,x), 



(32) 



where 



F{h,x) 



sin ^(h) — sin 



h 



X 



(33) 



1 — ^ 

and where < a; < 1 and —l<h<l. Similarly, Qn{s) in Eq. (|2^) can be defined in terms 
of a generating function Pi(/i, s) of the form 

' hR 



F,{h,s) = Y.Qn{s)h'' 



n=0 VI - /i^P^ 

where \hR\ < 1. Qn{s) is then given by 



sin ^{hR) — sin' 



1 - hRjl - ^ 



(34) 



(35) 



Tn{s) in Eq. (^) can also be defined in terms of a generating function F2{h, s), 

' hRs 



F2{h,s) = J2Z,{s)h'' 



1 



1 - hRsJl - ^ 



(36) 



where \hRs\ < 1. T„(s) is then given by 



(37) 



It follows from the previous results that Pn{s) can be expressed in terms of two unique 
elementary functions, Fi{h,s) in Eq. ( p4D and F2{h,s) in Eq. (|36D , such that 

^ {mT m = o,s) _ .-li, {ly F,{h = 0, .) 



Puis) 



(38) 



It is convenient to summarize the different expressions we have obtained for the PDF of 
a uniform n-dimensional sphere with radius R: 



1. Integral representation: 



Puis) 



1-1 n — 1 



hB{l + \,\)R^ 



n 



(39) 



2. Normalized incomplete beta function representation: 



Puis) = n- 



R 



n 1- 



n 1 1 
" ' 2 + 2'2/' 



(40) 



3. Incomplete beta function representation: 



PJs) = n- 



(- + - - 

V2 ^ 2' 2 



4. Odd-integer finite series expansion representation (n = odd): 
PJs) = nx 



n-l 

2 (_iY 



■ . V (¥)' 

i?" (n-l)!!^2z + l,!(ii^_, 



2i?y 



2i+l' 



5. Even-integer finite series expansion representation (n = even): 



Pn{s) =nx 



R" 



2 (n- 2i)\\ 

-^o^" [7715] --22 



n-2i + l 

,2\ 



(41) 



(42) 



\2Rj TT^ (n-2i + l)!! 



R 



2i-2-n 



(43) 
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6. InGnite series expansion representation: 

Pn{s) 







(2.-1)!! (^)! 


/ \2j+l" 

i-(ft) 


^ 2i+l 


(2i)!! 




2n 


(2 ~^ 2' 2) -^'^ 





(44) 



7. Generating function representation I: 



n\ V (9/1 7/1=0 



(0 



Puis) 



sin ^{hRs) — sin ^ 



/ifo- A / 1 



l-hRs^ I- 



2n-^V2 ^ 2' 2''-"' 



(45) 



8. Generating function representation 11: 



1 



Puis) 



sin ^{hR) — sin ^ 



hR-. 1 



2n-^V 2 ^ 2' 2''-'^ 



(46) 



9. Hypergeonietric function representation: 
2n 



Puis) 



R oF^ ( — , — : — : 1 

* 2 2 2 2 



- 2-ri 



2' 2 



2' 2' 4R2 



(47) 



where 



(48) 



and 

y{x) = 2Fi{a,b;c;x) 

abx a{a + l)b{b + 1) a{a + l){a + 2)b{b + l){b + 2) x^ 
c(c+l) 2[^ c(c+l)(c + 2) 

is one of the solutions for the hypergeometric equation 

a;(l — x)y"{x) + [c — {a + b + l)x] y\x) — aby{x) = 0. 



(49) 



(50) 



Using the previous results one can obtain a number of identities and recursion relations for 
the probability density functions Pn{s), as we discuss in Appendix 0. 
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III. SPHERICALLY SYMMETRIC DENSITY DISTRIBUTIONS 



In this section we generalize the previous results to the case of an n-dimensional sphere 
of radius R with a variable (but spherically symmetric) density distribution of the form 



p — p{r), where r = is measured from the center and 

xl + xl + --- + xl< R'^. 

As before we begin with the example of a circle (n = 2) and generalize to a sphere [n > 3) 
later. Following the derivation presented in the previous section, the positive x direction is 
chosen to specify the distribution of those s vectors that are aligned along the positive x 
direction. At this stage we must consider the differences between uniform and non-uniform 
density distributions. For a given s, if point 2 carries the density information p{x,y), then 
point 1 should have the density information p{x — s,y). It follows that to incorporate the 
effects of a spherically symmetric density distribution the following substitution should be 
made: 

p(f2) X p(fi) p{x, y) X p{x - s,y). (51) 

Since the density distributions considered are spherically symmetric, the probability of find- 
ing a given s in any orientation is still proportional to 27rs. The PDF P2{s) can then be 
expressed in the form 

p / N ^ M^) X Ms) + n2{s)) 

' /o^^n,(s) x(ni(s)+ 712(5)) ds' ^ ' 



where 



Udis) = 27rs, 



ni(s)=/ dx p{x,y) X p{x - s,y)dy, 

Js-R J—s/R^-{x-s\^ 



(x—s 



n2{s)= dx p{x,y) x p{x - s,y)dy. (53) 

•/f j-VR^-x^ 

Substituting x — s = x' and using p{x,y) = p{—x,y) it can be shown that ni(s) = n2(s). 
The expression for P2{s) can then be simplified to read 
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^ _ s dx /^^g^p(x, y) X p{x - s, y)dy ^^^^ 
lo^ (s dx I^^j^^ p{x, y) X p(x - s,y)dy^ ds 



The formalism leading to Eq. (54) can be extended to a 3-dimensional sphere of radius R. 
For a given s, the 2- axis is chosen arbitrarily as our reference axis to examine the distribution 
of those s vectors that are aligned along the positive z direction. If the density at point 2 is 
p{x, y, z), then the density at point 1 will be p{x, y,z — s). In analogy with the 2-dimensional 
case discussed above, the expression in Eq. (|5TD must be replaced by 



p(f2) X p(fi) = p{x, y, z) X p{x, y,z- s). (55) 

Since the density distributions considered here are spherically symmetric, the probability of 
finding a given s in any orientation is proportional to Atts^. Hence Psis) can be expressed 

as 



s'^l^dzj^^^^dxj^^^^^p{x,y,z)xp{x,y,z-s)dy ^^^^ 
J^^ (^s^ J^dzJ^^^^dxJ^^^^^p{x,y,z) X p{x,y,z- s)dy) ds 



J^dx3 j^^-^^^dXi J^^JZ^^dX2P (Xi, X2, X3) p (Xi, X2, X3) 



lo^ s'^ds J? dxs J^^l—^ dxi J^' ^^/'/'a dx2p (xi, Xa, X3) p (xi, X2, x'3) 



(57) 



where X3 = X3 — s. 

Up to this point our discussion has been completely general. To continue we next evaluate 
Psis) for a 3-dimensional sphere of radius R using two different spherically symmetric density 
distributions. Consider first 

P(r) = J^r\ (68) 

where N = An r'^p{r)dr, and r is measured from the center of the spherical distribution. 
Combining Eqs. (|57D and (|58|) we find 



^3(^) - - ^^R^~ 16 R^ + 448:R^^" ^ ^ 

A plot of Eq. ( p^ when R = 1, along with the corresponding Monte Carlo results, is shown 
in Fig. |. The second spherically symmetric distribution we consider is 

13 



p(r) 



N 



4vr i-f]/23 



(60) 



where N = Arc r'^p{r)dr, < a < 1, and r is measured from the center. Combining 
Eqs. (13) and (ID we find 

15(35 - 42a + 15a2)s2 225(1 - a)^^^^ 15as^ 



+ 



7(5 - 3a)2i?3 4(5 - 3a)2i?4 (5 _ Sa)R^ 

75(l + 6a-3a2)^5 ^5^^? 45aV 



(61) 



16(5 - 3a)2i?6 8(5 - 3ayR^ 448(5 - 3ayR^° ' 
A plot of Eq. (|6l|) when i? = a = 1 is shown in Fig. along with the corresponding Monte 
Carlo results. 

A general formula for the probability density function for an n-dimensional sphere with 
radius R having a spherically symmetric density distribution can be derived from the pre- 
vious results. We find 



2 _^_R2_^2 



/ 



yi2 ^2 ... I"' 2 



p (X) p (X') dXn- 



Jo 



2R 



2?!— 1 



Js dXn J 



dxi 



I 



:^p(X)p(X') dXn-l 



ds 



(62) 



where 



p(X) = p(xi,a;2,X3, 

p(X') = p(xi,X2,X3, 



(63) 



Another density distribution we wish to study is a Gaussian. As an example, consider 
the case of an n-dimensional sphere of radius R 00 with a Gaussian density distribution 
p(r) given by 

N 



where 



(27r)20-' 



-e ^ 



2 



(64) 



TT 2 

= lim n — 7 ^ , 

r (f + 1) -^0 



R 



Pn{^)'^^ dr. 



(65) 
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In Eq. (^5]) r is measured from the center of the spherical distribution and the integral is 
over all space. Recall that 

' x'^e-^dx = 2^r f ^) a"+^ (66) 



Combining Eqs. ( |5^ and the PDF for an n-dimensional sphere in an infinite space 
with a Gaussian density distribution can be expressed as 

P„(s) = lim = s'^-^e-^. (67) 



For n = 3, Eq. (|67|) agrees with the result obtained earlier in Ref. [T^. Finally, we note that 
the maximum probability, denoted by Smax, occurs at 



Smax = J2{n- 1) a. (68) 



IV. ARBITRARY DENSITY DISTRIBUTIONS 

We consider in this section the probability density functions for an n-dimensional sphere 
of radius R having an arbitrary density distribution, 

p = p(X) = p(xi,X2, ■ ■ ■ ,Xn), (69) 

where 

xj + xl + --- + xl<R\ (70) 

The proportionality factors, 27rs {n = 2) and Atts'^ {n = 3), cannot be applied here directly 
because the density function is not spherically symmetric. For a given s, each direction of s 
carries different information specified by the density distribution function. 

We begin with a circle of radius R and the conventional notation for polar coordinates, 
X = rcos(f) and y = rsin0. For a given s = r 2 — fi, the PDF, P2{s), is proportional 
to Ag X p(r2) X p{fi), where Ag is the overlapping area between the original circle and a 
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second identical one whose center is shifted to s, as described in the previous sections. In 
2-dimensional space s can be characterized by an angle (j) in the range < < 27r. One can 
understand the new features that arise for a non-uniform density distribution by referring 
back to Fig. |^. In the case of a uniform density distribution the picture formed by the vector 
s extending between Li and L2 is unchanged by a rotation of the entire pattern about the 
positive X-axis. However, for a non-uniform distribution the effect of such a rotation is to 
shift the vectors into a new region for which the density of points is not the same as it was 
initially. Stated another way, for a fixed Is*! the shape of the overlapping area or volume is 
the same, but will contain a different fraction of the points depending on the orientation of 
s. To deal with this effect, one can rotate the coordinate system so that the pattern remains 
as shown in Fig. |l], but with an appropriately transformed density distribution. To specify 
this transformation, we associate s with a rotation operator R(s ) such that the direction 
of s is the new x direction where s = x' , x' ■ x = cos0, x' ■ y = sin0, y' ■ x = — sin0 and 
y' ■ y = cos 0. We utilize the transformation matrix for a 2-dimensional rotation 



-^2x2(0) 



cos (p sm < 
— sin 6 cos I 



(71) 



to describe R(s). Notice that -R2x2(0) is an orthogonal matrix which satisfies -^2x2(0) = 
-^2x2(0)' where T denotes the transpose. Recall that the functional form of the density dis- 
tribution in Eq. ( 169|) is written in the original coordinate system. The inverse transformation 
matrix, -^2x2(0) = -^2x2(0) 5 should be used to transmit the correct density information to 
the new coordinate system s {x'). 

It is convenient to introduce the following general notations. 



X' = i?L2(0)X, 







x' 








_y'_ 





cos </> — sm I 
sin 6 cos d) 



(72) 
(73) 



and 



X" = i?J,2(0)(X-S) 



(74) 
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x" 




COS0 


— sin 




X — s 


y" . 




sin (j) 


COS0 




y 



(75) 



where 



X 



X 




x' 




x" 


, s = 


s 


y 




y' _ 




y" _ 








(76) 



We can then express the PDF for a circle of radius R with an arbitrary density distribution 
as 



s C # /f dx L^^g. p(X')p(X")dy 



/o k /r /f dx r^;^ p(X')p(X")dy 



r2^ j-R fVR^-x'- 

2 



(is 



(77) 



where 



p(X') = p(cos X — sin y , sin x + cos (p y) , 
p(X") = p(cos (x — s) — sin y, sin (x — s) + cos (f) y) . 



(78) 
(79) 



As an example , consider a circle of radius it! and non-uniform density distribution p{x, y) 
given by 



37ri?io 

where is a normalization constant 

rR 



. , 640iV 4 4 640iV 8 4^.4 
p{x,y) = T^—p^x y = - — zrn^r cos (psm 



(80) 



TV 



dx / p(x,y)dy. 



(81) 



For this example the 2-dimensional PDF P2{s) is then given by 

o9 



875 s 500 7400 400 

+ 



81 ^2 3 21 i?6 3 ^8 

-[/l(«)+/2(«)-/3(«)] 
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TT 



1750 s 1000 14800 800 s 



7 



+ 



+ 



+ 



+ 20 



(82) 



where 



17 



Ms) 



14875 92500 

162 ^ 
260315 5^° 



553985 



243 1701 R^' 



113693 



2509 



as 



10206 47g28 i?i6 142884 i?20' 

2725 1438825 s^^ ggigg gie 



+ 



+ 



1134 142884 285768 



(83) 
(84) 
(85) 



Figure]^ exhibits P2{s) when R = 1, and illustrates the agreement between the Monte Carlo 
simulation and the analytical result given above. 

The preceding discussion can be extended to a 3-dimensional sphere of radius R with 
an arbitrary density distribution p = p{x, y, z), where x = r sin 6 cos (p, y = r sin 6 sin 0, and 
z = rcosO are the usual 3-dimensional spherical coordinates, and x'^ + y'^ + z'^ < R^. For 
a given s = f2 — fi, the PDF ^3(5) is proportional to p(r2) x p(^i) x Vs, where Vg is the 
overlapping volume between the original sphere and a second identical one whose center is 
shifted to s. In 3-dimensional space s can be oriented at any angle between and 2tt, and 
the angle 9 can lie between and vr. A rotation matrix R^^^iO^cp) is used to represent the 
rotation operator R(s) associated with a given s such that 



i?3x3(^,0)=i?3x3(^)xi?3x3(0) 



COS 9 cos COS 9 sin — sin 9 

— sin cos 
sin 9 cos sin 6' sin cos 6* 



(86) 



where 



R3x3i9) 



i?3x3(0) 



COS — sin 9 

1 

sin 9 cos 

cos sin 

— sin cos 

1 



(87) 



(88) 



We observe the following: 



1. The rotation matrices R 



3x3 1 



)) -^8x3(6'), and -R3x3(0) are orthogonal so that 
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(89) 



2. The purpose of -R3x3((/>) is to transform the coordinate system from (x, |/, z) to a second 
coordinate system [x^^y^^z^) given by 



(90) 



3. The purpose of R3x3{d) is to transform the coordinate system from {x^,y^,z^) to a 
third coordinate system {x^,y^,z'^) where 



x^ 




COS0 


sin0 







X 


y' 




— sin (p 


COS0 







y 


z' 










1 




z 







x"^ 








y' 













cos 6 — sin ^ 

1 
sin 9 cos 9 



X 



Define the following notations, 



X!^Rl,,{9A)X 



(91) 



(92) 



where 



X" = i?J,3(e»(X-S) 

= i?J,3(0)i?J,3(^)(X-S), 





X 




x' 




x" 







x = 


y 




y' 




y" 


, s = 







z 




z' 




z" 




1 



The PDF ^3(5) can then be expressed in the form 

S In aiLiUUU In UUJ I s U,^ I 



s lo sin 9d9 d(p /. dz dx J 



(93) 



(94) 



gLp(X')p(X")dy 



/o^^ sin 9d9 C dct> J? dz dx 11^^ p(X')p(X")dy 



ds 



(95) 
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where 



x' = cos 9 cos (j)x — sin (py + sm9 cos 2;, 
y' = cos ^ sin (j)x + cos y + sin ^ sin (f) z, 

z' = — sm9 X + cos9 z, (96) 

pi^") = pix",y",z"), 

x" = cos 6* cos (p X — sin (j) y + sin 6* cos (j) {z — s), 
y" = cos ^ sin X + cos (py + sm9 sin (2; — s) , 
2;" = — sin 6* a; + cos 9 {z — s). (97) 

As an example, consider a 3- dimensional sphere of radius R and non-uniform density 
distribution p{x, y, z) given by 

945A^ 2 2 2 945A^ r . 4 ^ 2n 2,-2, /^on 
p{x,y,z) = y z = j:^r sm 6' cos 6' cos 0sm 0, (98) 

where N is the normalizing factor 



Then 

1701 _ 25515 8505 _ 8505 567 6237 / 

~ Tis"^ ^ 572 ^^l43'^~^08"^'^Tr^~lor^^ ^ 
201285 181629 16443 s^^ 6075 10899 

^ 9152 ~ 18304 R^^ ^ 6656 ~ 18304 + 585728 

Figure |^ is the plot of Psis) for i? = 1, and illustrates the agreement between Monte Carlo 
simulation and the analytical result. 

We can extend the discussion to a 4-dimensional sphere of radius R and arbitrary density 
distribution function, 

p = p(xi,X2,X3,X4), (101) 
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where 



(102) 



The 4-dimensional hyperspherical coordinates |jl3[ that are a generahzation of the conven- 
tional 3- dimensional spherical coordinates are defined as follows: 



Xi 


= r sin 62 


sin 61 cos 


X2 


= r sin 62 


sin 61 sin 


X3 


= r sin 62 


cos6'i, 


X4 


= r cos 02 


5 



and 



r = + X2 + xl + xl, 



61 = tan 



2^3 



02 = tan" 



X4 



_i X2 
tan — , 

2^1 



where 



(103) 



(104) 



(105) 



0<r <R, 0< 91,02 <n, 0<(p<2n, 

and the volume element dV is given by 

dV = dxidx 2(1x^(1x4^ = sin^ 62 sin9idrd92d9id(f). (106) 

The representation of the rotation operator R(s) for a given s is a 4-dimensional rotation 
matrix 



R. 



■4x4lC'2, 9i, 



-^4x4(^2) X -^4x4(^1) X -R4x4(0) 

COS 9i COS (f) COS 9i sin (j) 

— sin (j) cos (p 

cos 6*2 sin 6*1 cos cos 6^2 sin 6*1 sin 1 

sin 92 sin cos sin ^2 sin 9i sin < 



- sin 



cos 92 COS 9i — sin 6*2 

sin 6*2 cos 6*1 cos 6^2 



(107) 
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where 



-R4x4(6'2) 



-^4x4(^1) 



i?4x4(0) 



10 
10 

COS ^2 —sin ^2 
sin 6*2 cos ^2 

cos^i — sin^i 

10 

sin^i cos^i 

1 

COS0 sin0 

— sin0 COS0 

10 

1 



It is convenient to introduce the notations, 



X' — -RJx4(^2, Oi, 0)X, 

= -RIx4(0)-R4x4(^i)-R4x4(^2)X, 



X" = i?rx4(^2,^l,0)(X-S), 

= -RIx4(0)-R4x4(^i)-^4x4(^2)(X - 

where 



Xi 










X2 




x'. 




X2 


X3 










X4 




x^ 




■t'4 



The PDF ^4(5) can then be written as 
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p /X _ I [O2: 0l:4']x I [X1:X2, X3, X4] X p{X.') X p{X.") 

I^''{s^I[e2,e,,<P]xJ[x,,X2,X3,x4xp{X')Xp{X")}ds' ^ ' 



where 



[92,9i,(j)]= sm^92d92 sm9id9i d(p, (115) 

Jo Jo Jo 

J [-^'-^'-3,X4] ^ dx, dx, l^^^^^ dX2 dXs, (116) 

x[ — cos cos (pxi — sin + cos 92 sin cos + sin ^2 sin 9i cos 0a;4, 

^2 = cos 9i sin + cos + cos ^2 sin ^1 sin + sin ^2 sin ^1 sin (px^, 
x'^ — — sin ^i^i + cos ^2 cos ^i^s + sin 92 cos ^^1X4, 

X4 = — sin ^22^3 + cos 92X4, (117) 



p(X")=p(a;;',4',4',a;:), 

x'/ = cos 9i cos — sin 0X2 + cos ^2 sin 9i cos + sin ^2 sin ^1 cos 0(x4 — s) , 
X2 = cos 9i sin 0xi + cos 0X2 + cos 92 sin 9i sin 0X3 + sin 92 sin 6*1 sin 0(x4 — s) , 
X3 = — sin ^1X1 + cos 92 cos ^1X3 + sin ^2 cos 9i (x4 — s) , 

X4 = — sin 6'2X3 + cos 6'2(x4 — s). (118) 

As an example, consider a 4-dimensional sphere of radius R and non-uniform density 
distribution 

. 327V 4 12Ar 4 . 4„ . 4„ 4 , 
p(xi,X2,X3,X4) = ;^j^^Xi = -^jy^r sm 6'2sm e^icos 0, (119) 

where the normalizing factor N is given by 



N= dXi ^ dX2 / , ciX3 / ^ p(Xi,X2,X3,X4)dX4. (120) 

J-R i-^/R2Z^ J^./R2_^2^2 }_ /R2_^1 1 I 



-^R^'-xl J-^R^-x\-xl J-^R^-x\- 

Then 
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1/196 114 28 s*^ 4 5^0 2 s^'^ I s^^\ r— ^ 



TT 



3 i?6 ' 5 /28 ' 15/^10 5/^12 9^14 45^16 



1 / 112 



7r \ 3 

Figure | is the plot of Pa{s) for i? = 1, and illustrates the agreement between the Monte 
Carlo simulation and the analytical result. 

We turn next to the general case of an n-dimensional sphere of radius R and arbitrary 
density distribution, 



p = p(xi,X2, ■ ■ ■ (122) 



where 



xl + xl + ■■■ + xl<R^. (123) 
Define the following n-dimensional spherical coordinates Xi, . . . , x„ ||14||, 



Xi = r sin 9n-2 sin ^„_3 sin 62 sin 9i cos < 

X2 = r sin 6'„_2 sin 6'„_3 sin 6*2 sin 6^1 sin ^ 

= r sin On-2 sin 6'„_3 sin 6^2 cos 61, 



Xi = r sin On-2 sin 6'„_3 sin 6'j_i cos 



2, 



Xn_2 = r sin 6'„_2 sin 6'„_3 cos 6'„_4, 
x„_i = rsin6'„_2Cos6'„_3, 
x„ = rcos6'n_2, (124) 



where 



dV = dxidx2 ■ ■ ■ dxn 

= r"-^ sin"-2 9^-2 sin"-=^ 6^-3 sin^ 62 sin 61 drrf^„_2rf^n-3 rf^2rf^irf0, (125) 
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and 



0<r<R, 0< 91,02, ,0n-3,0n-2<Tr, 0<(/.<27r. (126) 
The rotation operator R(s) for a given s is 

(0). (127) 



The matrix Rnxni(p) appearing in Eq. (|127|) has the following elements: Rii{4>) = coscp. 



Ri2{4>) = sin0, -R2i(0) = -sin0 , R22{4>) = cos0 , and Rim{4>) = ^im where l,m ^ 1,2 
and \ < l,m < n. The matrix Rnxni^i) has the following elements: i?ii(^^i) = cos^i, 
-Ri3(^i) = -sin 6*1, i?3i(6'i) = sin 611, RssiOi) = cos6'i, and Rim{Oi) = ^im where l,m ^ 1,3 
and 1 < l,m < n. The matrix Rnxni^i) has the following elements: 1(6*4) = cos6'j, 

Ri+i^i+2{0i) = -sm9i, Ri+2,i+i{9i) = sinOi, Ri+2,i+2{di) = cos9i, and Rim{di) = ^im where 
2 <i <n — 2,l,m ^ i + l,i + 2, and 1 < l,m < n. Notice that RnxniOn-2, ^n-s, ■■ - 62,61, 0) 
has the following properties: 

1. Rnxn{0n-2, On-3, " " " 6*2, 6*1, (f)) is an orthogonal matrix such that Rnxn = ^nxn- 

2. The 1st row matrix elements are Ru = cos 6*1 cos 0, R12 = cos ^1 sin 0, Ri^ = — sin 6*1, 
and Rij = for 4 < j < n. 

3. The 2nd row matrix elements are R21 = — sin 0, R22 = cos 0, and i?2j = for 3 < j < 
n. 

4. The ith row matrix elements, where 3 < i < n — 1, are Ru = cos0j_i x xi[i], Ri2 = 
cos6'j_i X X2[i], Rim = cos6'j_i X Xm[i] for 1 < m < i, Ru = cos6'j_i x Xi[i], Ra+i = 
— sin^j_i, and Rij = for i + 2 < j < n, where Xm[i] is the mth component of the 
z-dimensional Cartesian coordinate system in the representation of the i-dimensional 
spherical coordinate system for a unit vector. Some examples are 0:3 [3] = cos^i, 
X3 [4] = sin 02 cos 9i , X3 [5] = sin 63 sin 62 cos 9i , and X3 [6] = sin 6*4 sin 63 sin 6*2 cos 9i . 
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5. The nth row matrix elements are Rnj — Xjln] for 1 < j < n, where Xj[n\ is the 
jth component of the n-dimensional Cartesian coordinate in the representation of the 
n-dimensional spherical coordinate system for a unit vector. Some examples are 



Xi[n] 


= sin 9n- 


_2 Sm 9n-3 ■ ■ ■ 


• • • sin 92 sin 9i cos 0, 




X2[n\ 


— sin 


-2 sin 9n-3 ■ ■ ■ 


• • - sin ^2 sin^i sin0, 




xsin] 


= sin 9n- 


_2 sm 9n-3 ■ ■ ■ 


• • - sin ^2 cos^i, 




Xn[n\ 


— cos 9n- 


-2. 




(128) 



The final master probability density function formula Pnis) for an n-dimensional sphere 
of radius R and arbitrary density distribution has the following mathematical representation: 

where 



j [9n-2, 9n-3: 92,9^^] = [ Siu^'^ ^n-2rf^n-2 [ Siu""^ 9^-3(^9^-3 

r sm^92d92 r sm9id9i d(t), (130) 

Jo Jo Jo 



\[X1,X2, ,Xn\= f dXn ^ dxi 



dX2 



-y/R^-xl J_^R2_^2_^2 

dXn-2 / dXn-l. (131) 

p(X')=p(a;;,4, x^), (132) 

p(X")=p«,4', 4). (133) 

Additionally, we introduce the following notations: 

= -Rnxn(^n-2, ^2, ^1, 0)X 

= -^nxn(^)-Rnxn(^l)-^nxn(^2) -R^xn(^n-3)-Rnxn(^n-2)X, (134) 
X" = -R^xn(^n-2, ^2 , ^1 , 0) (X - S) 

= -^nxn(^)-Rnxn(^l)-^nxn(^2) -R^xn(^n-3)-Rnxn(^n-2) (X - S), (135) 
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where 



X 



Xi 
X2 



Xr>. 



X' 



X2 



X" 



X. 



,11 

n-l 
II 



X 



n 



(136) 



The technique for generating random points within an n-dimensional sphere having an 
arbitrary density distribution will be discussed elsewhere [00] • 



V. APPLICATIONS 



A. mth moment 



We first calculate the mth moment (s"*) for the case of ra-dimensional uniform sphere, 
where 

I-2R 

(s^) = / s"'P^{s)ds, (137) 







and m is a positive integer. Evidently (s™) gives the expectation (average) value of the mth 
power of the distance between two independent random points generated inside a uniform 
n-dimensional sphere. By utilizing the function C (a; m, n) defined in Eq. (0), we can write 

(s"^) as 



C(2i?;0,n) \m + nj 5 (f + |, i 

where we have used the following identity: 

C{a; m, n) = —B - - + - W?" - - - + - IR" 



2m + n V2 2 2/ 2m + ra [jk) \2 2 2 

1 (27?)™+" /n 1 m n 1\ ^ , , 

+ -- — 5. ,^2 - + - + —,- + - LR". 139 

2m + n ^ V2 2 2'2 2y ^ ' 
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Furthermore we can rewrite Eq. ( p.38| ) in terms of the gamma function by replacing the beta 
functions as follows: 

r(f + f + i)r(n + i) 



{s"") = {2Ry 



n 



m 



n 



+ ^^r(n + i + f)r(f + i) 



2 r(n + m + i)r(|) 



+ n) r(f + f)r(n+i + f; 

The results in Eqs. ( |140| ) and ( |141|) are identical to those given in Refs. 
extended to evaluate and we find 

on-m B (- + - ^ + i — 



(140) 



1411 



1 

am 



n-m RfB + ii 

V2 2' 2 



They can be 



(142) 



where 

m < n — 1. 

Combining Eqs. (|138D and (|142|) , the mth moment (s™) has the general form 



(143) 



n 



n + m 



13 ( n \ I. B_|_i_|_Hi 
-°V2~'~2'2~'~2~'~ 2 



R ( — -\- - - 

V 2 2' 2 



(144) 



where 



m = -in- 1), -(n - 2) , . . . , -2, -1, 0, 1, 2, 



(145) 



Following is a short list of (s"^) in 3 dimensions: (l/s^)=9/4, (l/s)=6/5, (s)=36/35, 
(s2)=6/5, (s3)=32/21, (s^)=72/35, and (s5)=32/ll, where the radius R has been set to 
unity. 

Additionally, (s™) can be evaluated for a sphere having a Gaussian density distribution 
and radius R ^ oo. From Eq. (^) we have, 



r2R r f ^) 

= hm / s"^P^is)ds = {2ar^^. 
R^oo Jo r I - 1 



(146) 

In some applications involving low-energy interactions among nucleons the lower limit 
(zero) should be replaced by the hard-core radius Vc = 0.5 x 10^^^ cm |TH]. In such cases 
the expressions for P„(s) and (s™) assume the form: 
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/f dx - ^ C'(2i?; 0, n) - C(r,; 0, n) 



(147) 



and 



is ; J^^ S I:^n{S}as ^^^^ , 



(148) 



where 



H {R, Tc] m, n) 



{2RY+^ r /n In 1 m\ 

r U ^ 2' 2 + 2 + 2"y' ~ ^(ft) 

r"+™ \^fl n l\ ^ /I n 1 

5 -,- + - -S,, ,2 -,- + - 

n + m [ V2' 2 27 ift) V2' 2 2 



n In 1 mV 
^'2 + 2'2 + 2 + Yj. 



(149) 



and m is an integer. 



B. Coulomb Self-Energy of a Collection of Charges 

As an application of the preceding formahsm we evaluate the electrostatic energy Wn 
of a collection of Z charges in n dimensions by applying geometric probability techniques. 
Consider first the familiar case of a spherical charge distribution in 3 dimensions. For each 
pair of charges the potential energy due to the Coulomb interaction in Gaussian units is 

v(\r2-n\)^j^^, (150) 

1^2 - ri| 

where Cq is the elementary charge {el/hc = 1/137). Hence if we assume that the charges in 
each pair are uniformly distributed within the same spherical volume of radius R then the 
average Coulomb energy C/3 of each pair of charges is 

Us^e' !Lp,{s)ds^^^. (151) 
Jo s o K 

For a collection of Z charges there are Z{Z — l)/2 such pairs, and hence the total Coulomb 
energy W3 is 

W^-^^^Us^lz{Z-lfj. (152) 
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For n > 3 the Coulomb potential energy between two charges has the general form |T6 



Vn{\r2-r,\)= (153) 
r2 - ri 



where g„ is a suitably defined charge with appropriate dimensions. Hence Eq. ( |151| ) gener- 
alizes to 

r2R I / 1 \ 

Ur^ = ^nl ^2Pn{s)ds = qli-^). (154) 

Using the results of Eq. (|144|) with m = n — 2 we then find 

;^„.Z,Z_l)^illMl^.fj!L)£(^. ,155) 

^ \2 ^ 21 2 J 

For n very large Wn assumes the limiting form 



limW^n = 2 ^ (156) 



If the density distribution p„(r) of the charges is Gaussian rather than uniform, where 

Pn{r) = T^r^e-^S , (157) 
(27r) 2 an 



then 



1 o2 

Wn = T^^. (158) 



In the limit n —>■ oo, 



Wn = -e"/2n-"/2^. (159) 



These results are of interest in the context of recent work on modifications to the Newto- 
nian inverse-square law arising from the existence of extra spatial dimensions . It is well 
known that the gravitational interaction is weaker in a space with n > 3 spatial dimensions, 
this weakness being manifested by an inverse-power-law for the force F„, F„ oc l/r""^. 
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C. Neutrino-Pair Exchange Interactions 



A second example of interest is the z/i/-exchange (neutrino-pair exchange) contribution 
jlSl^Il to the self energy of a nucleus or a neutron star. For two point masses the 2-body 
potential energy is given by 



Vup{\fi-fj\) 



l5 ' 



(160) 



4vr^ |f 

where and aj are coupling constants which characterize the strength of the neutrino 
coupling to fermions i and j {i, j = electron, proton, or neutron). In the standard model |^ 

1 



+ 2 sin^ Ow 



2 sin^ e 



w 



0.964 
0.036 



1 

an - 

In contrast to the Coulomb interaction, the functional form of V^p(r) cannot be determined 
in a space of arbitrary dimensions on the basis of a general argument utilizing Gauss' law. 
Hence we restrict our attention here to 3 spatial dimensions and consider the case of a sphere 
of radius R containing neutrons. For the case of a uniform density distribution we then 
find 



N{N -I) 



9 



+ 



9 



3r^ 



Gl ^ ?,N{N - l)Gl 



2rm^ Ar,R^ 8R^ 16R^ / An^ 



;i61) 



where Vc is the hard-core radius. The analogous result for a Gaussian density distribution is 



r (0, rll^a^) 



N{N -l)G 



32(T%7/2 



(162) 



J.2 4^2 

where r(a, h) is the incomplete gamma function. The expression in Eq. (|161|) agrees with 
the result obtained in Refs. pTlESi. 



D. Neutron Star Models 

Another application of current interest is the self-energy of neutron star arising from the 
exchange of vv pairs. Here we evaluate the probability density functions in 3 dimensions 
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for neutron stars with a multiple- shell uniform density distribution, which is what is typi- 
cally assumed in neutron star models For illustrative purposes, we discuss spherically 
symmetric models with spherical 2, 3, and 4 shells, where for simplicity we assume shells 
of equal thickness. Some other multiple-shell models and their ra-dimensional probability 
density functions can be found in 

For a 2-shell model with a uniform density in each shell define p = pi for < r < R/2 
and p = p2 for R/2 < r < R, where pi and p2 are constants and r is measured from the 
center of the neutron star in 3 dimensions. Using the preceding formalism we can show that 
the PDF has 4 different functional forms specified by 4 regions: 

1. < s < ii? : 

^ ^ 24(p^ + 7pl)s' _ 36{pj - 2pip2 + 5pl)s^ Ujpj - 2p,p2 + 2pl)s' 

(Pi + 7p2)' R' (Pi + lp2f R' (Pi + lp2f R' ' ^ ' 

2. \R<s<R: 



81(pi-p2)p2S , 24pis2 36pi(pi + 3p2)s3 

-T SlSl — — — ,0 + 



2 (pi + 7p2)' (p, + 7p2) + ^p^Y ^4 

I2pls^ 



(Pl + 7p2)'i?6' 



3. R<s<lR: 



81 (pi - P2)P2S , 24(9pi - P2)P2S^ 36(5pi - P2)P2S' 



2 (pi + 7p2)' R^ (pi + 7p2)' R^ (pi + 7p2)' R" 
12{2p,-p2)p2s' 



(pi + 7p2)'i?6 



4. |i? < s < 2i? : 



(164) 



(165) 



192pig^ 144pi53 , 12pls' 

Psis) = r5 r5 9 . 166 

(Pi + 7p2)' R' (pi + 7p2)' R^ (pi + 7p2)' ^ ^ 



We observe that the PDFs defined in adjacent regions are continuous across the boundaries 
separating the regions. 
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For a 3-shell model of a sphere of radius in 3 dimensions, with a different uniform 
density in each shell, define p = pi for < r < i?/3, p = p2 for -R/3 < r < 2i?/3, and 
p = P3 for 2i?/3 < r < R, where pi, p2, and pa are constants and r is measured from the 
center of the neutron star. In this case the PDF has 6 different functional forms specified 
by 6 regions: 

1. 0<s<lR: 

81(p2 + 7p2 + I9p2)s2 729(p2 - 2p^p, + 5pl - Sp^Ps + I3pl)s^ 



(Pl + 7p2 + 19p3)2i?3 4(pi + 7p2 + 19p3)2i?4 

2187(p^ - 2pip2 + 2pi - 2p2P3 + 2p^)g5 

16(pi + 7p2 + 19p3)2i?6 



(167) 



2. ii? < s < fit: : 



81(p2-p3)(9pl -9p2 + 25p3)g 81(pf + 7pip2 - 7piP3 + 26p2P3)g' 

"^^'^ 8(pi + 7p2 + 19p3)2/?2 (p, + 7p2 + 19p3)2i?3 

729(p2 + 3pip2 - 5pip3 + 10p2P3)s3 2187(p? - 2pip3 + 2p2P3)s'^' ,^ 

7^7 TTF. — xopfi , l-LoSj 



4(pi + 7p2 + 19p3)2i?4 16(pi + 7p2 + 19p3)2i?6 



3. |i? < s < : 



p . , 81(9pip2 - 9p^ + 55pip3 - 30p2P3 - 25pi)5 81(9pip2 - p^ + 19pip3)g" 

"^^'^ 8(pi + 7p2 + 19p3)2i?2 (p, + 7p2 + 19p3)2i?3 

729(5pip2 - p2 + 5pip3)s3 2187(2pi-p2)p25^ ,^ 



4(pi + 7p2 + 19p3)2i?4 16(pi + 7p2 + 19p3)2i?6' 



4. < s < |i? : 



. , 81(64pi-39p2-25p3)p3g 81(8pi + 28piP3 - 9p2P3)g' 

'^^^ 8(pi + 7p2 + 19p3)2i?2 (pi + 7p2 + 19p3)2i?3 

729(4p2 + 10pip3 - 5p2P3)s3 2187(pi + 2pip3 - 2p2P3)/ 



4(pi + 7p2 + 19p3)2i?4 16(pi + 7p2 + 19p3)2i?6 



5. |i? < s < |i? : 



2025(p2 - P3)P3S , 81(35p2-8p3)p3s' 

^3[S) - 



8(pi + 7p2 + 19p3)2i?2 (p, + 7p2 + 19p3)2i?3 

729(13p2 - 4p3)p3S^ ^ 2187(2p2 - P3)p3g^ ^^^^^ 



4(pi + 7p2 + 19p3)2i?4 16(pi + 7p2 + 19p3)2i?6' 
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6. lR<s<2R: 

2187p^s2 6561p^s3 2187p^s5 



(Pl + + 19P3)'^=' 4(pi + 7p2 + 19p3)2i?4 16(pi + 7p2 + 19p3)2i?6 ' 

(172) 

As in the previous case, the various functional forms for ^3(5) are continuous across the 
boundaries separating the regions. 

The 3-shell model is of interest since actual models of neutron stars often invoke a 3- 
shell picture 0J^. We note to start with that although our 3-shell model assumes that all 
shells have the same thickness (shell radii of r, 2r, and 3r), we can relax this assumption 
by substituting r — > ri, 2r ^ r2, and 3r — > r3, where ri, r2, and are arbitrary numbers. 
Similarly, the densities pi, p2, and p3 can also assume arbitrary values, so that the results 
of Eqs. ( |167| )-( PT2| ) can be applied to any realistic neutron star model [0. 

The preceding formalism can be extended to any number of shells. For a 4-shell model 
define p = pi for < r < R/A, p = p2 for R/4 <r<R/2,p = p3 for R/2 < r < 3R/A, and 
p = Pi for 3-R/4 < r < R. Here pi, p2, ps, and p4 are constants, and r is measured from the 
center of the neutron star. The PDF has 8 different functional forms specified by 8 regions: 



1. R<s<lR: 



192ipl + 7pl + 19pl + 37pl)s' 

(pi + 7p2 + 19p3 + 37p4)2i?3 

576(pf - 2pip2 + 5pi - 8p2P3 + 13p^ - 18p3P4 + 25pl)s^ 

(pi + 7p2 + 19p3 + 37p4)2i?4 
768(p^ - 2pip2 + 2pl - 2p2P3 + 2pi - 2p3P4 + 2pl)s^ 

(pi + 7p2 + 19p3 + 37p4)2/?6 ' ^ ) 

2. \R<s<\R: 

p( ^ 18(9pip2 - 9p^ - 9pip3 + 34p2P3 - 25p^ - 25p2P4 + 74p3P4 - 49p|)g 

(pi + 7p2 + 19p3 + 37p4)2i?2 

192(p^ + 7pip2 - 7piP3 + 26p2P3 - 19p2P4 + 56p3P4)g^ 
(pi + 7p2 + 19p3 + 37p4)2i?3 
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576{pI + 3pip2 - 5piP3 + IOP2P3 - 13p2P4 + 2 0p3p4)g^ 

(174) 



(pi + 7p2 + 19p3 + 37p4)'i?^ 
768 (p^ - 2piP3 + 2p2P3 - 2p2P4 + 2p3P4)s^ 



(Pl + 7p2 + 19P3 + 37p4)2i?6 



3. lR<s<lR: 

18(9p^ - 9pip2 - 55pip3 + 30p2P3 + 25p^)g 
'^^^ (pi + 7p2 + 19p3 + 37p4)2i?2 

18((31pipi - 18.3p2Pi + 70p3Pi + 19p.^J.s 

^ (Pl + 7p2 + 19p3 + 37p4)2i?2 

192(9pip2 - pj + 19pip3 - 26pip4 + 63p2P4)g' 

(pi + 7p2 + 19p3 + 37p4)'i?-^ 
576(5pip2 - P2 + 5pip3 - 10piP4 + 17p2P4)g^ 

(Pi + 7p2 + 19p3 + 37p4)2^^ 
768(2pip3 - P2 ^ 2pip4 + 2p2P4)5^ 
(pi + 7p2 + 19p3 + 37p4)2i?6 ' 



4. lR<s<R: 

18((31pip3 - 39p2P3 - 25pf^ + IGlpiPi - 12p2Pi - 70p3pi - lOp^g 

'^^'^ ~ (Pl + 7p2 + 19p3 + 37p4)2i?2 

192(8pi + 28pip3 - 9p2P3 + 37pip4)g^ 
^ (Pi + 7p2 + 19p3 + 37p4)'^^ 
576(4p| + 10pip3 - 5p2P3 + 7piP4)g'^ 

(Pl + 7p2 + 19p3 + 37p4)2i?4 

768(p^ + 2pip3 - 2p2P3)g' 

(Pl + 7p2 + 19p3 + 37p4)2i?6' 



18(25p2P3 - 25p^ + 225pipi - lOGpaPi - 70p3Pi - 19pl)g 

(Pl + 7p2 + 19p3 + 37p4)2i?2 

192(35p2P3 - 8pl + 65pip4 - 28p2P4)g^ 

(pi + 7p2 + 19p3 + 37p^i)^R^ 
576(13p2P3 - 4pi + 17pip4 - 10p2P4)g^ 

(pi + 7p2 + 19p3 + 37p4)2i?4 
768(2p2P3 - P3 + 2piP4 - 2p2p4)s^ 

(Pl + 7p2 + 19p3 + 37p4)2i?6 ' 
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6. f i? < s < |i? 



'^^^ (Pl + 7p2 + 19p3 + 37p4)2i?2 (p, + 7p2 + 19p3 + 37p4)2i?3 

576(9pi + 20p2P4 - 13p3P4)s3 768(pi + 2p2P4 - 2p3P4)s5 



(pi + 7p2 + 19p3 + 37p4)2i?4 (p, + + I9p3 + 37p4)2/?6 ^ 



7. |i? < s < : 



o , , 882(p3 - p4)p4S , 192(91p3 - 27p4)p4S 



(pi + 7p2 + 19p3 + 37p4)2i?2 (p, + 7p2 + 19p3 + 37p4)2i?3 
576(25p3 - 9p4)p4S^ 768(2p2 - P4)P4S^ 



+ 7--^77^^''^w::''^'^^ (179) 



(pi + 7p2 + 19p3 + 37p4)2i?4 (p, + 7p2 + 19p3 + 37p4)2i?6 



8. \R<s<2R: 



12288p2s2 9216pls3 



(pi + 7p2 + 19p3 + 37p4)2i?3 (p, + 7p2 + 19p3 + 37p4)2i?4 
768pls^ 

^(Pl + 7p2 + 19p3 + 37p4)2i?6- 



(180) 



A statistically tested random number generator and a special simulation technique have 
been employed to numerically simulate the probability density functions, and the numerical 
results agree with our analytical results as discussed in Refs RH]. 



E. Geometric Probability Constants 

We can also apply the preceding geometric probability techniques to evaluate the expec- 
tation values of fi2 ■ 7^23 in n dimensions, where fi2 = ^2 — ri, r23 = rs — 'r2, and ri, r2 and r^, 
are three independent points produced randomly inside an n-dimensional uniform sphere of 
radius R. The quantity (ri2 ■ 'r23)„ is one of the geometric probability constants of interest 
and has important applications in physics 0,0. In 3 dimensions. 



(ri2 ■ r23)3 = ((^2 - xi) (X3 - X2) + (?/2 - yi) {ys - 2/2) + {Z2 - Zi) {zs - Z2)) . (181) 
Following Refs [^i. 
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where 



J% X V (n - 1, - x"^) dx n 



V{n,R) = —, ^. (183) 

r(f + i) 



We can verify Eq. (182) by applying the geometric probabihty techniques directly so that 



1 /"SH n 
{ri2 ■ r23)„ = "2 I s'P4s)ds = '^R'- (184) 



Following is a short list of {fi2 




for selected values of n: 




{ri2 ■ 


^23)1 = 






(ri2 ■ 


^23)2 = 


2 ' 




{ri2 ■ 


^23)3 = 








^23)4 = 






{ri2 ■ 






{ri2 ■ 


^23)5 = 




Notice that as n becomes larj 


le (fi2 


■^23)^ - 


If we 



(185) 
(186) 
(187) 
(188) 
(189) 



dimensional spherical space with a Gaussian density distribution as defined in Eq. (|6^) and 
radius — > 00, then 

(ri2 ■ r23)„ = -na\ (190) 

An interesting application of our formalism in biomedical physics is the calculation of the 
dynamic probability density function P(s, t) in 3 dimensions, as discussed in Ref. [§. 



VI. CONCLUSIONS 

A new formalism has been presented in this paper for evaluating the analytical probabil- 
ity density function of distance between two randomly sampled points in an n-dimensional 
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sphere with an arbitrary density. We illustrate the power of this new geometric probabil- 
ity technique by demonstrating how n-dimensional integrals can be reduced 1-dimensional 
integrals, even in the presence of Vc- We have shown that the classical result of Eq. ( |28D 
can be rederived from a simple intuitive picture, and that it can also be extended to an 
n-dimensional sphere with a non-uniform density distribution. Several density examples 
(Eqs. (^8|), (|60|) , (|80|), (p8|), and ( |119| )) were presented, and the analytical results were 
shown to be in agreement with Monte Carlo simulations. Applications to a variety of phys- 
ical systems, such as neutron star models, have also been discussed. 
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APPENDIX A: GEOMETRY OF INTERSECTING CIRCLES 

Here we briefly show why circles Oi and O2 are identical with the center of Oi located 
at (0, 0) and O2 located at (s, 0) as shown in Fig. |[ Following the discussion in Sec. [II A| , 
define a Cartesian coordinate system for points 1, 2, 3, and 4 as 



= (^s/2, ^R^~sy4j, (Al) 

{X2,y2) = -^R^ - syAj , (A2) 

(2^3,1/3) =(s-i?, 0), (A3) 

(x4,i/4) ={R,0). (A4) 



Assume that the equation for the circle O2 is 
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I n2 I 2 2 

[x — a) + y = r , 



(A5) 



where (a,0) is the center and r is the radius. Inserting Eqs. ( |A1| ), (^), and ( [A3[ ) into 
Eq. ( [A5D which expresses the fact that the circle O2 contains points 1, 2, and 3. If s 7^ 2i?, 
then the only solution is a = s and r = R. This result means that circles Oi and O2 have 
identical radii and the center of O2 is located at (s, 0). 

APPENDIX B: RECURSION RELATIONS OF P^is) 

We present in this Appendix some recursion relations for the probability density functions 
Pn(s) which follow from the results in the text. 



PL(s) = ^^PJs) 



n 



n-l 

.2 \ — 



1 



(Bl) 



n — 3 

2 \ — 



45 (f + i, I) V 4i?2 



(B2) 



\-^\p:is)-'^i2-^-^^PLis^ 



n — 



/2n - 1\ s 



\ A J 



Pn{s) = 0. (B3) 



- 1) (A) 



0. 



(B4) 



n + l 

2 \ — 



n R^ ' Tx (n + l)!!i?"+3 \^ 4P2^ 



n = even 



(B5) 



n s 



n + 2 s 
1 n!! 



n 



2P2 



Pn-2(S) 



2 



1 + ra 



4P2 



4P2 



even) . 



(B6) 



Pn+2{S) 



n + 2 



■PJs 



1 (n + 2)!! s"+2 



.2 \ 2 



n P2 2(n + l)!!i?"+3 I 4P2 



(n = odd). 



(B7) 
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n 



n+2 



2(n + l)!!i?"+i 



n s 
1 + n 



Pn-2{S) 



„ , , n + 2 ^ , . 



r,n+2 



n 



2Pn{s) 



n 

n + 2~s^^ 
1 



n+2 



n s 
n-2W 



n-2 



1 + n- 



n 
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FIGURES 
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(a) (b) 

FIG. 2. (a) Locus of points in a circle of radius R separated by a vector s = sx. For each point 
fi in Ai, there is a unique point r2 in ^2 such that s = r2 — H- ^2 is the overlap between the circle 
and the region obtained by translating the line L2 given by x = s/2 to the +x direction. Ai can 
be obtained by shifting A2 a distance s along the +x direction, (b) A^ is the overlap between the 
circle and the region obtained by translating the line Li given by x = — s/2 along the —x direction. 
A4 can be obtained by shifting ^3 a distance s along the x direction. 
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s 




FIG. 3. Geometric interpretation of the probability density function P2{s). From Figs. |T] and 
the probabihty that two points are separated by a distance s {\s\ = \sx\) in a circle of radius R 
is proportional to the total shaded area, IJ IJ ^4 U ^2 1 in the circle Cq. As can be seen from 
this figure, this shaded area is given by the overlap of Co with two identical circles Ci and C2 as 
shown, where Co is the circle + = i?^, Ci is (x — s)^ + = R^, and C2 is (x + s)^ + = R^. 
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analytical 


o 


simulation 



1.2 - 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

S 

FIG. 5. PDF for a 3-dimensional sphere of radius R = 1 and p oc 1 — r^. 
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FIG. 9. A simple diagram showing circle Oi and O2 are identical with the center of Oi located 
at (0,0) and O2 at (s,0). 
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